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■»>.^ 1 Einstein (1936) predicted the phenomenon presently known as gravitational lensing (GL). A prime 

• feature of GL is the magnification, because of the gravitational field, of the star visible surface as seen 

^L(' from a distant observer. We show here that nonlinear electrodynamics (NLED) modifies in a funda- 

mental basis Einstein's general relativistic (GR) original derivation. The effect becomes apparent by 
studying the light propagation from a strongly magnetic (B) pulsar (SMP). Unlike its GR counterpart, 
the photon dynamics in NLED leads to a new effective GL, which depends also on the B-field per- 
^ ' meating the pulsar. The apparent radius of a SMP appears then unexpectedly diminished, by a large 

factor, as compared to the classical Einstein's prediction. This may prove very crucial in determining 
physical properties of high B-field stars from their X-ray emission. 
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1. INTRODUCTION 

Einstein's general theory of relativity (GTR) has proved to be one of the most success- 
ful physical theories ever formulated. As a theory of the gravitational interaction, it has 
a number of predictions definitely corroborated by experiments or astronomical observa- 
tions, which makes it the correct gravity theory here-to-fore. After readdressing an earlier 
study, Einstein (1936)^came up with the prediction of the gravitational lensing effect. The 
idea is that the gravitational field produced by a massive astrophysical object, the Sun for 
instance, can act as a convergent lens able to deviate light rays flying-by. Observations of 
far away quasars and galaxies and more familiar total solar eclipses confirm the reality of 
this phenomenon. More than half of the star's surface may be seen by a distant observer 
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This means that a photon emitted at a given colatitude on the star's surface, reaching the ob- 
server at infinity, must be emitted at a smaller angle with respect to the normal at that point. 
Because of this effect, a distant compact star must then appear to any observer larger then 
it actually is. The relation between the apparent radius i?oo of a spherical (non-rotating) 
star, as seen by a distant observer, and its physical radius R is 

This relation is obtained from the photon trajectory in polar coordinates (r, 9) by using 
a Schwarzschild metric 

V r 

Here Rs = 2GAd /(? is the Schwarzschild radius of the star of mass M. As is clear 
from this line-element no effects from physical fields other than the gravitational have 
been taken into account in prescribing the star apparent size. This paper shows that the 
phenomenon changes in a fundamental fashion if NLED is called into play. The analysis 
below, based on well-known nonUnear Lagrangeans, as the Heisenberg-Euler (H-E,^ and 
the exact Born & Infeld (1934) (B-lEl), proves that when NLED is included to describe the 
photon dynamics, the trajectory depends on the background B-field pervading the star. We 
stress from the very beginning that the results below are obtained upon an idealization of 
the i?-field structure. In a real star, the B-field is neither purely radial nor constant, so that 
the apparent radius hardly goes to zero at high B. 

We stress in passing that this theory is not a new theory of gravitation. In our formalism 
Einstein's general relativity remains definitively unaltered. Only the photons propagating 
out of the ultramagnetized neutron stars "see" the effective metric that appears once one 
incorporates NLED into a self-consistent physical description of light propagation from 
those stars. That NLED is a theory akeady tested in a number of experiments is demon- 
strated by the laboratory experiment described by Burke, Field, Horton-Smith, et al. (1 997) 
1^ as quoted below. Our theory attempts to call to the attention of the astrophysics com- 
munity the fundamental changes that should appear in determining physical properties as 
radius, mass, and more crucial: the equation of state, of neutron star pulsars endowed with 
extremely strong magnetic fields, since a property as the star radius (which can be esti- 
mated from the x-ray variability of very short period of the star emission (light curve), or 
by measuring the gravitational redshift of absorption or emission lines received from the 
star surface, or etc.) can say too much about the stellar structure and its equation of state. 
A fundamental property that encodes basic information of the physical constituents of the 
neutron star matter 

2. Apparent radius in NLED 

According to quantum electrodynamics (QED: see Delphenich 200 3 Bl for a complete re- 
view on NLED and QED) a vacuum has nonlinear properties 12151 which affect the pho- 
ton propagation. A noticeable advance in the realization of this theoretical prediction has 
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been provided by Burke, Field, Horton-Smith, et al. (1997)^, who demonstrated ex- 
perimentally that the inelastic scattering of laser photons by gamma-rays is definitely 
a nonlinear phenomenon. The propagation of photons in NEED has been examined by 
several authors (Bialynicka-Birula & Bialynicka-Birula 1970; Salazar et al. 1989; Di- 
etrich & Gies 1998; De Lorenci 2000; Denisov & Svertilov 2003; Alarcon 1981; Gar- 
cia & Plebanski 1989; Boillat 1970; Vazquez et al. 1969; Torrence 1984; Garcia 1984) 
| 7 | 8 | 9 | 10 | 17 | 1 1 | 12 | 14 | 15 | 16 EI in the geometric optics approximation, it was shown by 
Novello et al. (2000)1^, and Novello & Salim (2001) that when the photon propa- 
gation is identified with the propagation of discontinuities of the EM field in a nonlinear 
regime, a remarkable feature affloats: the discontinuities propagate along null geodesies of 
an effective geometry which depends on the EM field on the background. This means that 
the NEED interaction can be geometrized, in analogy to gravity in GTR. A key outcome 
of this formalism is introduced in this paper. 



2.1. Euler-Heisenberg approach 
The H-E'2lLagrangean 

L{F^G)^-\f+^^^{f^ + \g'^ , (3) 

9 2 

where = is a quantum parameter, is a gauge invariance description of the 

photon propagation in NEED that uses two invariants 

F = F,,,F^'', G^F;,F^'' = K^^,''"Fa.pF^'', (4) 

constructed upon the Maxwell EM tensor and its dual ^. The method of characteristic 
surfaces or shock waves followed here (introduced by Hadamard (1903)'—') can be applied 
to any field theory having hyperbolic field equations, including electrodynamics. Hence, 
the result in Ea.(l50> (see appendix), implies two possible paths of propagation or polariza- 
tion modes, according to the double solution il±. Using Eq.Q to compute the Lagrangean 
derivatives Lf,Lff, Lqg, one arrives to the couple of effective metrics (first order in /i) 

<?,^ff% = 9^"" + ^tiF^^F^- , g^^_ = g^'' + UtiF^^F^^ . (5) 

In terms of electric {E) and magnetic {E) fields the tensor can be written as 

F^,, = E^V, - E,V^, + ii^^u'^'^V^Bp, (6) 

where is the normalized {V^^V^ — 1, Eq.Q) velocity of the reference frame where 
the fields are measured. Our main concern here is the behavior of photons, in the NEED 
context, emitted from a highly magnetized neutron star. Corotating charges in the pulsar 
magnetosphere or a rotating magnetic dipole lead to induced i?-fields in the star surface. 
We consider here slowly rotating neutron stars in order that the i?-field contribution could 



^The attentive reader must notice that this first order approximation is valid only for B-fields smaller than Bq = 
™M = 4.41 X 10^^ G (Schwinger's critical B-field). 
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be neglected. Since E"' — 0, the above expression simplifies to: i^^^ = Vi-tiy^^^aBp. Its 
self-product reads 

pt^ap^^ ^ _ ^2 ^^M^ _ V'^V) , (7) 

where ~ B'^B^. After discarding a nonphysical conformal factor the effective 
metrics become 

ll^ = 5'^'' - . I b^^B'^ + l/'^l-- . (8) 



' 1-/3252 



The inverse (covariant) metric of Eq.(|S) is obtained from the relation ^^jyg"" — 5^°" , 
which then reads 

gf.^ = .9m- + PlB'^B'' - PlB^V'V . (9) 

By we mean = lAfi and /3i = 8/i. In order to pursue our calculations we as- 
sume a radial magnetic field. This is clearly not realistic, but is an approximate description 
of the field geometry near the polar caps of a magnetized neutron star Thus 

B^^rn^, l^ = V^r5l, V^^^oSl, (10) 
and consequently one arrives to the effective metric components 

5So± = (1 - /3±S') 500 , gfr^ = {l-~PlB^)grr ■ (11) 

Notice the structural similarity of both metric components. From Eq. dl It is straightfor- 
ward to rewrite the expression for the apparent radius of the spherical star in Eq.Q in the 
putative background metric given by Eq.(|2j. After doing so, one can compute the ratio of 
the star area with H-E NEED effects included to the area without (Einstein's derivation), 
as seen from the distant observer This quantity is represented here by M. ^ 



2.1.1. Qualitative description of the visible ( apparent) star surface 

As stated above, "A/"" is intended to be a quantity describing the ratio between the area of 
the star as seen by an observer at infinity that receives the actual traveling luminous ray, 
and the actual physical area of the neutron star, both projected on the plane of the sky. This 
ratio of areas can be consistently figured out by the following qualitative reasoning. 
In a general relativistic framework, the local neutron star luminosity is defined as 

= 4^i^2^sB^,^ (12) 



''We stress here that the NLED effects depend on the effective geometry, which in turn depends upon the strength 
of the magnetic field. Although the magnetic field _B is a solution of the NLED field equations, to compute 
the ratio "A/"" is only needed to know the local (on the star surface) magnitude or modulus of the field. In the 
particular case of the computation of the null geodesies discussed below, is indeed necessary to know the analytic 
expression for the field as a function of the coordinates, in particular of the radial r-coordinate. Unfortunately, 
such a solution for the magnetic field structure (or dependence on the coordinates) is only known for the linear 
(Maxwell's) case. 
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Fig. 1. H-E NLED to GR area ratio, as seen from a distant observer, for an arbitrary M/ R ratio. Notice that 
the ratio grows down to zero near the QED critical field _B ~ 4.4 X lO^^ G. This means that no visible area is 
seen from the SMP. This "anomalous" behavior invalidates the Denisov et al. results; which advogate for crucial 
changes in using the H-E approximation (see discussion below). 

where Ri, is the neutron star radius and is its temperature. Both measured locally 
at the star surface. The luminosity measured by the observer at infinity, which carries the 
information on the spacetime metric around the neutron star, is defined 

Lao^L^[l-^^, (13) 
where ^1 — is the tt-metric component of a Schwarzschild spacetime , i.e., 

.,.= (l-f). (.4) 

Notice that grr — j^- On the other hand, the luminosity at infinity can be written as 

Loo = inRl^asBirn^. (15) 

Thus, the temperature measured by an observer at infinity relates to the actual star 
temperature via 



- (^1 - j . (16) 

Therefore, after manipulating all these relations one arrives to the relation between both 
the star radii (local and at infinity) 

i?oo - i?4.9rr)'/'. (17) 

In an equivalent approach, once one realizes that the star luminosity is modified by the 
spacetime geometry induced by the matter distribution, one can readily verify that in the 
case of NLED its geometric effects on the photon propagation from the ultra magnetized 
neutron star leads to an expression that is similar to the one above. Hence, one can write. 
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taking into account the general relativistic effects superposed to the NLED effects, the local 
neutron star luminosity as given by 

= AiTB?asiiTt (18) 

The luminosity measured by the observer at infinity, which carries the information on 
the spacetime metric as modified by NLED, as in the Heisenberg-Euler (H-E) NLED, now 
reads 

Rs 



1 - PiB' 



or equivalentiy. 



where the term 



L^J,f=4n{RtlffasB{Tn^ 



(19) 



(20) 



1 



now substitutes the tt-metric component of 
the pure Schwarzschild spacetime presented above. Notice that in this case the effective 
tt-metric component reads 

Rs 

This impUes that 



1 



1 - (3lB^ 



(21) 



g-^f = {grr)\l-~0iB 



1 - /32 B2 



9tt 



Thus, the temperature measured by an observer at infinity now reads 



R* 



1/2 



1/2 



Once again, after manipulating these relations one obtains 

R'" = R* 



(22) 



(23) 



(24) 



From these relations one can write the area ratio "J\f" that estimates the magnitude of 
the change in the apparent visible area of the star as seen from infinity compared to the 
actual star area, as measured by a local observer : 

Areaiocai 47r(i?*)2 



AreaNLED 4^(^^/)2- 
Hence, by substituting the corresponding terms in these last equations one arrives to 



(25) 



geff 
yrr 



1 - plB^ 



(26) 



As one can notice, in this new expression for the area ratio there is no any functional 
of the r coordinate but rather only a function of the magnetic field strength at the surface 
of the neutron star. This is the reason of why we took B as a constant in the numerical 
calculation leading to the plots in Figures 1, and 2. 
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2.2. Born-Infeld NLED 

The propagation of light from hypermagnetized neutron stars can also be viewed within 
the framework of the Born-lnfeld Lagrangean 

i = ( W 1 + ^ - G2 - 1 I (27) 




where h"^ — -h = jm^c? — > = 9.8 x 10^^ e.s.u. As is well known, this is 
an exception al Lagr angean. One of its remarkable properties is that it does not exhibit 

IT1 "1 o' 

birefringence ° . In this case, the deduction we present in the appendix fails since the 
quantities fi^, with i — 1,2,3; vanish identically. We need hence to go back to Eg. ( 1461471 : 
and substituting in these equations the form for L{F^ G) of Ea.( l27t to obtain a unique 
characteristic equation, which then leads to the effective metric 

5c"ff + (28) 

Using the self-product of the tensor F^i, given by 0, and noting that in our case 
F = FujyF^'^ = 2B^, the effective metric then reads (see reference^^for details) 



"62 



By computing the inverse metric via g^^gf,^ = S^^, the covariant form of this effective 
metric is obtained as 

2B^ /62 2_B2 /h'^ 

= 5m- - (252/^,2 + 1) + (252/52 + 1) ^'^''^ • (^^^ 

After recalling the relations for and given in Ea.(llO>. one can verify that the 
covariant rr effective metric component is then written as 

1 



i/rr .yi 



(252/52 + 1) 



1 + 252/52 



Qrr ■ (31) 



By following steps similar to the H-E case, the area ratio is presented in Fig. 2. Although 
our hypothesis about the field geometry maximizes the NLED effect on the observed prop- 
erties of the star, they should properly be taken into account when analyzing the X-ray 
emission. (This point will be considered in a forthcoming paper). It should also be rele- 
vant in pondering the effects on supernova dynamics of photons radiated either by highly 
magnetized proto-neutron stars (Miralles et al. (2002) or stellar-mass black holes en- 
shrouded by super strong 5-fields (see van Putten et al. (2004) 26). 



3. Photon trajectories as seen from a distance 

A self-consistent treatment of this problem does require the study of the light propaga- 
tion in the neutron star spacetime geometry modified by NLED. In NLED the Lagrangean 
defined by Eq.(|2j becomes 
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Fig. 2. H-E NLED to GR area ratio, as seen from a distant observer, for an arbitrary M/ R ratio. Notice that the 
ratio grows down to zero near the QED critical field B ^ 4A X lO^^ G. This means that no visible area is seen 
from the SMP. 



Rs\ f dtV n{B) f dr^"^ 



r 



dXj (1- ^) V*^^ 



where ij,{B) = -TTprhuw- resulting Euler-Lagrange equations read (with li" 



d 



— (-2g,f w*) = u* = Eo, (33) 



d 

dX 

d_ 

dX 



— {2g'eeu') = geeu' - Lo, (34) 
{2gllu^) = gllu^ = h, (35) 



where h is defined as the impact parameter of the photon propagating from the NS 
surface to an observer at r — > oo, and Lo a constant, which may be put equal to zero. 

As the photons travel along null geodesies, the Euler-Lagrange equation for the rr- 
component of the metric can be obtained from the (first integral) relation 

gtfiu'f+gf^u^f+geeiu'f+gf^u^f ^ . (36) 

In analogy with the planetary system one can chose a particular "orbital" plane at = 
7r/2, so that = and sin^ 0=1, and thus g^^ = r^. By making use of the standard 
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change of variables {u ~ 1/r,' — d/dr, and defining A = 2i3^/6^) one can write the 
photon propagation equation as 

^""(aTt) (^)/iV = 0, (37) 

which one can in turn write as {Eq, h constants) 

[ur + = § + Rsu' + f 2§ -u^- RsuA A+^A^. (38) 

After performing the derivative of equation (I38> with respect to the angular coordinate 
(j), and taking into account the fact that the solution u = constant can be discarted by 
an analogous (physical) reason as it is done in the case of a pure photon propagation in a 
Schwarzschild geometry, one arrives to 

It" + (1 + A)u = ^Rs (1 -A)u^ + 

This equation is highly nonlinear and quite hard to solve by analytical means. (Notice 
that Ea. (l39> is quite similar to that one obtains in the case of the photon propagation in 
a Schwarzschild spacetime. Nonetheless, we stress that in order to find such a solution; 
an analogous procedure as the one usually pursued in general relativity {the perturbation 
approach) does not work in the present case, since for the extremely large magnetic field 
we are considering the induced NEED effects are of the same order of magnitude as those 
purely obtained from gravitation.) Besides, one needs to know the function B{r, </)), which 
is a solution of the Born-lnfeld equations for this problem. Perhaps a numerical solution 
can be found, but that would require extra work that we think would not add too much to 
the new results we are presenting in this version of the paper. (A full detailed solution will 
be given elsewhere). 

Eas.( l38l3'9t shows very clearly that NEED does modify the standard propagation of 
photons as compared to that one in a pure Schwarzschild gravitational fieldl^^l as described 
in most textbooks of GTR. A first look at Ea. (l38> suggests that very large B-fields do 
reduce the effective star's visible area. 



1 + RsU 



2^(1 + A)u-' 



dA 

dr 



^- (39) 



4. Discussion and conclusion 

The nonlinear electrodynamics is a well stablished theory of the electromagnetic interac- 
tion, which overcomes the failure of Maxwell's theory in describing the stability of atomic 

o 

structures. Since its introduction in physics by Euler and Heisenberg it has gained the 
status of a theory with a strong experimental support (for a recent complete review see 
Delphenich (2003)ii!). Despite being known since long, its very key features have become 
recognized only in the late years (see references above). It had been shown in the last few 
years (Novello et al2000; Novello & Salim 2001) that the force acting along the photons 
trajectory can be geometrized in such a way that in an effective metric g'^(,-^ = g^^i, +.9^^^^^ 
the photons follow geodesic paths in that effective geometry (Novello et al2000; Novello & 
Salim 2001) '^^'^^ ' . In particular, the effects of nonlinear electrodynamics in the physics of 
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Strongly ma gnetized neutron stars have been recently studied by Mosquera Cuesta & Salim 
(2003; 2004)1221241 It was shown there that for extremely supercritical magnetic fields non- 
linear electrodynamics effects force photons to propagate along accelerated curves, in such 
a way that the surface gravitational redshift of emission lines from hypermagnetized neu- 
tron stars, as measured by a distant observer, is significantly modified. Such a new effect 
turns out to be of the same order of magnitude of the one produced by the pure gravitational 
field. 

In this paper, we present a new astrophysical application of the peculiar effects brought 
by NLED, which is in timatel y related to the one previously discussed by (Mosquera Cuesta 
& Salim 2003; 2004) ' -^ l '. We show that the famous gravitational lensing effect, origi- 
nally introduced by Einstein, is also modified when the electromagnetic (nonlinear) effects 
are taken into the depiction of the physics around strongly magnetized neutron stars. Such a 
dynamics is described by using both the Euler-Heisenberg and the Born-Infeld approaches 
to NLED. Incidentally, we show, by comparing Eqs. Jl 1> and ( I30t . that Eq. Jl U clearly ex- 
hibits a divergence. Such an unphysical behavior, in turn, invalidates the results presented 
by Denisov et al. in Ref.^^ obtained within the H-E approach, since they extended to the 
case of magnetars (neutron stars having B-fields ~ 10^^ G) their formula to estimate the 
light-ray bending angle beyond the QED S-field limit. That is inconsistent. Besides, both 
approaches: the H-E and the B-1 NLED induce critical changes in the rr metric compo- 
nent that modifies the quantification of the lensing effect. More relevant, yet, the area of 
a putative star may appear to a distant observer nonphysically diminished in H-E NLED, 
whereas it is physically largely reduced, and may even "disappear", for fields B > 10^^ G 
in the B-1 approach. The impact of this on the star flux r nav be dramatic . Thence, the effect 
introduced by (Mosquera Cuesta & Salim 2003; 2004) I21I24I and this new here crucially 
alters the dynamics of photons from very high S-field stars, while entangles (makes it 
difficult) the inference of their physical properties. 

To the very end, what kind of observations can be performed in a search for this pe- 
culiar effect? As light is the prime messenger from the stars, one can look for variability 
of very short period in the x-ray (or 7-ray) light curves from most hypermagnetized pul- 
sars. A comparative analysis of the characteristics of those light curves with the ones from 
canonical pulsars may render the clues on the presence of NLED effects. A similar analysis 
can be pursued by using the measurements of the surface gravitational redshift of that kind 
of hypermagnetic stars. 

5. APPENDIX: The method of effective geometry 

Following Hadamard (1903)123 the surface of discontinuity ^ of the EM field is denoted by 
S. The field is continuous when crossing E, while its first derivative presents a finite dis- 
continuity. These properties are specified as follows: [-Fjii/lj^ = , [-F)ji/|a] ^ ~ //ji^^a , 
where the symbol — lim5^o+ ( J|s+(5 ~ J\^~s) represents the discontinuity of the 



Of course, the entire discussion onwards could alternatively be rephrased using concepts more familiar to the 
astronomy community as that of light rays used for describing propagation of EM waves in geometric optics. 
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arbitrary function J through the surface S. The tensor /^j, is called the discontinuity of the 
field, k\ — d\Yl is the propagation vector, and the symbol "|" stands for partial derivative. 

Here-after we investigate the effects of nonlinearities of very strong _B-fields in the 
evolution of EM waves; described onwards as the surface of discontinuity of the EM field 
(represented here-to-fore by Ffj,,^). Extremizing the Lagrangian with respect to the poten- 
tials Afj^ yields the following field equation: 

V^iLFFf"" + LgF*'"') = 0. (40) 

Besides this, we have the cyclic identity: 

V,F*^'' = ^ ^^.|" + ^"Mk + ^'^"IM = ■ (41) 
The field equation can be written explicitly as: 

LfV^F^'' + 2N''^"^>^V^Fa,f3 = , (42) 
where the tensor N is defined as 



j^^^af} ^ Lpj,F'"'F°'f^ + Lfg{F'"'F*"'^ (43) 
_|_ p*tiiypap-^ _j_ l^^^p*"'^ p*i^'^ _ 

Taking the discontinuities of the field equation we get: 



Uxk^ = -^N^'^'Upk^. (44) 
The discontinuity of the Bianchi identity yields: 

fapk^ + I^akp + fp-fka = 0. (45) 
From these equations we obtain (see^^lfor details) 

Xfc' = ■^F^-F\k,kr{LFFX + LgfX*) 
Lf 

- ^k\LFGX + LGGX*) (46) 
Lf 



^*fc2 ^ ^ptJ-^pr^j^^j^^^p^^^ ^ LgGX*) - 

Lf Lf 



2 



2F 

[LffX + LfgX*) + -j—k^iLFGX + LggX*) , (47) 
Lf 

where we introduce the notation: x = F""^ Jap, x* = F^^f"-^, fc^ = g^^'^k^k,^. As the 
H-E QED Lagrangean is not a functional of the product F * G, then one obtains 



Xk^ = -^F^-'F^k^kriLFFX) " ^k\+LGGX*) (48) 
Lf Lf 



I 
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= -^F'^''F\,Kkr{LGGX*) - ^k\LFFx) 
Lip Lp 

2F 

+ k\LGGX*) ■ (49) 
Lp 

We seek thus for a master relation representing the propagation of field discontinuities, 
which should be independent of the quantities fap, that is, independent of x and x* ■ There 
is a simple way to achieve such a goal. We firstly isolate the common term F^'^F'^^fc^fc^ 
which appears in both Eas.( l48l4'9t . Then, by assuming that k"^ ^ 0, the difference of 
these equations can be put in the form of an algebraic linear relation between x ™d x* 
given as: Vlix*'^ + ^2XX* + ^3X^ = 0, where we define: Vti = -^Lqq, and fl2 = 
Lgg + -^LppLggi ™d 5I3 — —^L^Fp- Solving the quadratic equation for x* we 
obtain: x* = ^^±X' with: fl± — ^'^2±y/n^---inin3 ^ ^^^^^ ^^^^ solution in Eqs.(l48l and 
( l49i . and assuming x 7^ 0' after some algebra one gets the dispersion relation 

3'"' - 4 , /^^^ r F^^F'^>) Kk. = • (50) 

Lp + Li\1±Lgg J 



Thence, one concludes that the discontinuities will follow geodesies in this effective 
metric gf^. 
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